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Abstract: We argue that it may be possible to reheat the universe after inflation driven | ^ 

by D-brane annihilation, due to the coupling of massless fields to the time-dependent 
tachyon condensate which describes the annihilation process. This mechanism can work 
if the original branes annihilate to a stable brane containing the standard model. Given 
reasonable assumptions about the shape of the tachyon background configuration and the 
size of the relevant extra dimension, the reheating can be efficient enough to overcome 
the problem of the universe being perpetually dominated by cold dark tachyon matter. In 
particular, reheating is most efficient when the final brane codimension is large, and when 
the derivatives of the tachyon background are large. 
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1. Introduction 

Because of the difficulty of finding direct laboratory probes of string theory, it is interest- 
ing to look for possible evidence from cosmology. Most notably, inflation may be tested 
more sensitively in the near future by the MAP ||^ and PLANCK [Q] observations of the 
cosmic microwave background radiation. Recently there has been significant progress in 
constructing stringy inflation models which make use of naturally occuring potentials be- 
tween D-branes to provide the false vacuum energy |^, |5|. However, these attempts have 
not adequately addressed the question of how reheating occurs after inflation. In fact, 
there are reasons to fear that reheating may be generically difficult to achieve in D-brane 
inflation. For this reason, we aim to propose a generic mechanism for reheating in such 
models, which is qualitatively different from reheating in ordinary field-theory inflation 
models, and which has the hope of being fairly robust. It is based on particle produc- 
tion in a time-varying background, which will occur even if the background motion is not 
oscillatory. 

Let us begin by describing the difficulty with reheating. In the simplest version of D- 
brane inflation, a parallel brane and antibrane begin with some separation between them 
in one of the extra dimensions required by string theory. Although parallel branes are 
supersymmetric and have no force between them, the brane-antibrane system breaks su- 
persymmetry so that there is an attractive force and hence a nonvanishing potential energy. 



It is the latter which drives inflation. Once the branes have reached a critical separation, 
they become unstable to annihilation. The instability is described by condensation of a 
tachyonic mode |^. Its low energy effective description is a field theory of a peculiar kind, 
whose Lagrangian has the form §1^ 

C = -Te~^"'-'F[-{d^Tf] (1.1) 

where 

Fix) - (1 2) 

^^''>~ 2T{2x) ' ^ ^ ^ 

which is determined for the superstring from the boundary string field theory (BSFT). Here 
T is the sum of the original brane tensions and a is of order the string length, a ~ ^foi' (the 
precise value to be used in the model we adopt will be discussed in section IV). The tachyon 
starts from the unstable maximum T = and rolls to T ^ oo. This process requires an 
infinite amount of time, during which the tachyon fluid has an equation of state identical 
to that of pressureless dust as T — > 1 ^] . 

(There have been numerous recent attempts to make use of the tachyon fluid for 
cosmology [pT| , either as the inflaton or as quintessence. Although these ideas might work 
if one had the freedom to change the form of the tachyon potential, the action which arises 
from string theory is not suitable for either purpose. Since the late-time equation of state 
is p = 0, the string theory tachyon does not provide accelerated expansion in the recent 
universe. At early times, its potential is too steep to satisfy the requirements of inflation. 



Constraints on tachyon cosmology and its shortcomings have been discussed in |12| and 

il.) 

Returning to our description of the endpoint of D-brane inflation, the situation is 



similar to that of hybrid inflation |14|, where the tachyon plays the role of the unstable 
direction in field space which allows for inflation to quickly end. The important difference 
is that in a normal hybrid inflation model, T would have a minimum at some finite value, 
e.g., due to a potential like A(|Tp — a^)^, and the oscillations of T around its minimum 
could give rise to reheating in the usual way, or the more efficient tachyonic preheating 



1 15]. But with the exponential potential in (|1.1D there can be no such oscillations. It 
thus appears that the universe will become immediately dominated by the cold tachyon 
fluid, and never resemble the big bang [l^. It is important to convert very nearly 100% 
of the energy stored in tachyon matter into radiation because the former redshifts more 
slowly than the latter. For example if we assume that the universe started initially with a 
temperature of at least 1 TeV for the purposes of baryogenesis, while matter domination 
begins at a temperature near 0.1 eV, then the ratio of energy density in tachyon matter 
versus that in radiation must have been no more than 10"^'^. 

To emphasize the difficulty of getting such efficient reheating in the present situation, 
let us contrast the rolling tachyon with an inflaton which is decaying through its oscillations 
into massless fermions via a coupling gcpipip. The probability of particle production in the 
background (p = 4>o sin(Mt) is proportional to the square of the transition matrix element 



^we use the metric signature (1,— 1,...,— 1) 



/ dtsm{Mt)e^''^* ~ 6{M - 2a;), where uj is the energy of the fermion or antifermion. The 
square of the delta function is understood in the usual fashion to be 6{M — 2uj) times the 
total time; in other words, oscillations lead to a constant rate of particle production. By 
simply waiting long enough, the energy stored in the oscillations will naturally be reduced 
to an exponentially small level. In the case of the rolling tachyon however, (j) is replaced by 
T which has the asymptotic behavior T ^ and the whole action is multiplied by a factor 
g-^/a or e~'^^/°'^ . The corresponding transition matrix element J dtte~*^°'^'^'^^^ is finite, 
with no delta function. As a result, it gives not a constant rate of particle production, 
but rather a finite number density of produced particles. It is not obvious that the energy 
stored in the tachyon fluid can be sufficiently reduced. 

Nevertheless, it is known that the tachyon couples to massless gauge fields; one form 



that has been suggested for the low energy theory is the Dirac-Born-Infeld action |16, 17| 



2 



C = -Te-^'/'^yi det(5^, + F^, - d^Td,T)\ (1.3) 

Because of its time dependence, we expect that some radiation will be produced by the 
rolling of the tachyon. It then becomes a quantitative question: can this effect be efficient 
enough to strongly deplete the energy density of the tachyon fluid, so that the universe 
starts out being dominated by radiation rather than cold dark matter? There is one 
immediate problem with this idea, however; the fact that the entire action is multiplied 
by the factor e"^^/'*^ means that the massless particles which are produced will not act 
like ordinary radiation |]l8|. Recent work has shown that these excitations have the same 
equation of state as the tachyon itself [^]. This is related to the fact that the system is 
annihilating to the closed string vacuum, which does not support any open string states 
like the gauge fields. 

An obvious way to circumvent the above difficulty is to assume that the branes anni- 
hilate not to the vacuum, but rather to branes of lower dimension [pO|]. This is a natural 
possibility since any initial fluctuations of the tachyon field where it changes sign will lead to 
the production of topological defects at the points where T = 0, via the Kibble mechanism 
Q. In the effective description, the lower dimensional branes are represented by solitonic 
configurations of the tachyon field. A D{p — 2) brane is a vortex of the complex tachyon 
field, whereas a D{p — 1) brane is a kink |0, |8|, ^ . These possibilities are described 
by the descent relations of Sen |Q. The stable descendant branes are able to support 
open string excitations of gauge fields, including those of the standard model. They will 
no longer decouple as a consequence of the rolling tachyon because the topological defect 
which prevents pins T = at the origin. 

In this paper we will set up a simplified model of particle production by tachyon 
condensation, which we hope captures the essential features of a more realistic Lagangian. 



^This form is not the same as what one derives from BSFT. However it has the right quahtative fea- 
tures, besides being simpler to work with. Ignoring for a moment the gauge field, it can be shown that 
\/7d£t^^)^r~^^"9^rrft7f^ = — for time-dependent configurations, so that the action vanishes as T — > 1. 



The same is true for the BSFT-derived function F(— T ) in (1.2). Moreover both functions behave like 



V (dxT)^ for spatially dependent profiles in the limit that dxT oo, whose relevance will be explained in 
the next section. 



The computational method developed here should carry over straightforwardly to more 
complicated situations. In section II we motivate an ansatz for the space- and time- 
dependent tachyon background which describes condensation to a brane of one dimension 
lower than the initial configuration. In section III we present the solutions for a gauge field 
in this classical tachyon background. Section IV describes how we compute the spectrum 
of particles produced during the early phase of the tachyon motion, and presents numerical 
results. We conclude in section V with the interpretation of these results and a discussion 
of how the calculation can be improved in future work. 

2. Tachyon background 

There arc two kinds of tachyonic solutions which have been described in the literature: 
(1) static solutions which are topological defects and represent lower dimensional branes, 
and (2) dynamical solutions which are spatially homogeneous and describe a cosmological 
fluid with vanishing pressure at late times. The mechanism we have in mind combines 
these two pictures by supposing that at t = the tachyon starts from (or very near to) the 
unstable equilibrium configuration T = throughout some number d of extra dimensions 
which will be transverse to the final descendant brane. We denote the spatial coordinates of 
these extra dimensions by .xi, . . . , x^- Starting at t = 0, the tachyon starts rolling towards 
T ^ oo for \xi\ sufficiently far from the center, whereas it is pinned to T = at = 0. By 
charge conservation, there must be at least one other position where T = also. The pair 
of defects represents a brane and an antibrane. For simplicity we will focus our attention 
on the half of the internal space near Xi = which contains the brane, since by symmetry 
the physics near the antibrane will be identical, for the purposes of particle production. 
At any given time, this half of the extra dimensional space will contain an interior region 
where the defect resides, and an exterior one where the tachyon is still rolling. We might 
expect the tachyon profile to evolve with time similarly to fig. 1(a). 

t 




X x=-L x=L 

Figure 1 . (a) Possible time evolution of the tachyon field T during condensation, (b) Our ansatz for 
the tachyon evolution in time and in the half of the extra dimensional space containing the defect 



representing the final state brane: region III contains a linear kink, II contains the homogeneous 
rolling field, and I is the unstable vacuum configuration. 

To find the real behavior of T(t, Xi) we should solve the equation of motion for T which 
follows from eq. (^]^). This is already a difficult numerical task in itself, and its results 
would surely complicate our next step, which will be to solve the gauge field equation 
of motion in the background T(t,Xi). Therefore we are going to satisfy ourselves with a 
simplified ansatz for the background and defer more detailed investigations to future work. 
Namely, we will consider a linear profile for the defect, and a purely time-dependent one 
in the exterior: 

T{x) = qx, \x\ < t/q (region III) (2.1) 

T{t,x) = aln(cosh(t/a))sgn(a;), \x\ >t/q (region II) (2-2) 
T = 0, t<0 (region I) (2.3) 

We have specialized to the case of a single codimension for the descendant brane {p = 1) 
and this dimension has been compactified by identifying the points x = 2L and x = — 2L, 
where an antikink is assumed to reside. We show only the half of the space —L < x < L 
which surrounds the kink. A parameter q describing the steepness of the tachyon kink 
profile has been introduced. The spacetime regions I-III are illustrated in figure 1(b). 



One may wonder to what extent the ansatz (2J-2^) reproduces a more realistic tachyon 
background. Let us first consider the rolling tachyon region (II). Eq. (2^) is in fact an exact 
solution not for the action (^]^) but rather for a mutilated version in which the argument 
of the exponential is linear: 



C = -Te-l^l/y I det((7^, - d^Td,T)\ (2.4) 

This closely resembles the BSFT result for the effective action from the bosonic string, 
except for having \T\ rather than T in the exponent. The bosonic string has an unphysical 
instability as T ^ — oo which we artificially remove by taking the absolute value. The 
behavior of the solutions to (2.4) is similar to that coming from a potential with a quadratic 



argument, which is the appropriate one for the superstring. For example, the late-time 
behavior corresponding to the Lagrangian £ = Te~"^^/'^^(1 — T^)^ is 

T{t) - t + ^!il_^(26)V(i"f)e-*V(a^(i-f)) + o(g-2tV(a^(i-f))) (2.5) 

assuming that b < 1. (This can most easily be derived by computing the corresponding 
Hamiltonian density and using energy conservation to get a first integral of the motion.) 



Similarly, if we expand (2.2) for large times, we get 

T{t) ^ t + oe-2*/« + Oie-^'/") (2.6) 

which resembles (^^). Moreover if we compute the relevant factor \/l — in the two 
cases (now taking 6 = 1/2), 

\/l-t2 = e"*'/'^' +0(e"2*V'^')^ e"*/"^ + 0(e-2Va) (2.7) 



respectively. In other words \/l — is just equal to the other prefactor e~*^/"^ or e~*/", 
as the case may be. This can easily be understood from the fact that the conserved 
Hamiltonian is V{T){1 — r^)-i/2^ hence V{T) and \/l — T^ must be proportional to each 
other. It will turn out that the simplified form e~*/'* makes it easier to solve for the gauge 
field in the next section, which is one of our motivations for adopting the action (p.4|). 

Next we discuss the ansatz for the kink profile, T = qx, in region III. This static 
profile is not a solution to the equations of motion^ except in the limit that q ^ oo. If 
we consider a generalized action of the form C = V{T){1 — T''^)^, then the linear profile is 
a solution only if = 1/(25 — 1). The limit g — > oo has been discussed in ||^, ^, where 
it was noted that the descendant brane resulting from tachyon condensation has the right 
tension to agree with string theory when g — > oo. In this limit the descendant brane looks 
like a genuinely lower dimensional object, whereas it would have a nonvanishing thickness 
(revealed by the energy density of the profile) for finite values of q. However a different 



method, that of level-truncation |22, leads to tachyon defect profiles which do have a 
nonzero width. In any case, there is no reason to believe that the kink thickness goes to 
zero immediately; it seems quite reasonable to assume that it will have some nonzero value 
initially, and possibly tend toward zero only as t ^ oo. We will assume that q remains 
approximately constant for < t < qL since this is the time during which particles are 
produced. If g ^ oo at later times, it will not significantly change our conclusions. 

A further issue is the shape of the spacetime boundary separating the static kink from 
the time-dependent solution. We have taken it to be linear, t = q\x\, but this implies that T 
is not even continuous at the interface except for t ^ a. It might seem more reasonable to 
deform the boundary such that e'^'^'l/'^ = e*/'^+e~*''". We could do so, but it would needlessly 
complicate the ansatz without even solving the problem of the tachyon Lagrangian being 
discontinuous at the interface, because the factor y^det(^^^i7— 7)^ )| = v^l - (a^r)2 is still 
not smooth: (d^T)'^ changes sign as well as magnitude across the interface: 



1 — e region II 

region III 



This is an indication that our ansatz (2J.-2!^) is defective; the true solution should have 



a smooth Lagrangian density at the transition. Nevertheless, it would greatly complicate 
the solution for the gauge fields in the next section if we tried to smooth out this behavior. 
We will instead compensate for the difficulties which arise from this oversimplification in 
another way, as will be explained. 

We found it simpler to consider condensation to a kink instead of a vortex configuration 
of the tachyon. The latter would be more realistic because of the fact that, assuming the 
parent branes are supersymmetric BPS states, a brane of only one dimension less is not, and 



^an exact solution can be obtained from that of region II by analytically continuing t — j- ix, giving 

' cos((I/ — x)/ay 



T = a sgn(x) In 



cos(L/a) 



where we have imposed continuity of T' at x = ±L and assumed that L/a < -k 12. The solutions for the 
gauge field are complicated in this background, so we prefer to use the linear kink profile for this paper. 



consequently not stable, whereas a brane whose dimensionality differs by an even number 
is BPS. To reduce the dimensionality by two, the tachyon should be a complex field which 
condenses to a codimension-two defect, a vortex. Although we originally wanted to treat 
this case, it is not clear how to write the action of the complex tachyon in a way which 
extends to the time-dependent configurations of region II. BSFT calculations show that for 
a winding configuration of the form 



the Lagrangian factorizes as 



T = qiXi+iq2X2 (2-10) 



C (X -F{ql)F{ql) (2.11) 



where F{x) = ^2r(2xj • "^^^ immediately obvious how to rewrite ( 2.11| ) in a Lorentz 
and gauge invariant way which would allow us to deduce the equations of motion for 
time-dependent homogeneous configuations. Explicit constructions involve the use of in- 
dependent tensors like d^T*d,^T and d^TdyT in matching powers of derivatives, so that a 
compact expression like ( |2.4| ) does not seem to emerge. We leave the consideration of these 
complications for future work. 

That being said, the kink configuration may still be physically relevant because of 
the possibility of descending from the original Dp-Dp pair in two steps, with the unstable 
D{p— 1) brane being a resonance through which the system passes on its way to the stable 



D{p — 2) endpoint [24 1. The transition from D{p — 1) to D{p — 2) would be described by 
the formation of a kink. 

Finally, let us consider the fact that the rolling phase of the tachyon field ends in a 
finite time within our ansatz. Fig. 1(b) shows that at late times we have eliminated the 
homogeneous condensate by fiat since at t = qL the entire bulk has been replaced by the 
static kink. We don't know whether this is the actual behavior, or if the homogeneous 
region persists, which could be the case if q grows with time: 

T{t,x) = q{t)\x\, \x\<t/q{t) (region III) (2.12) 

r(t,rc) = aln(cosh(t/a))sgn(x), \x\ > t/q{t) (region II) (2.13) 

If q{t) grows with t faster than linearly, then region II survives and the tachyon fiuid 
coexists with the final state brane at arbitrarily late times. We believe that the present 
calculation could give a reasonable approximation to the efficiency of particle production 
even in this case. The compactification length L will be replaced by the size of region III 
at the characteristic time scale when the fast roll phase ends {i.e. when T = 1 in the bulk). 
On the other hand, if q{t)/t ^ as t — > oo, then as long as the bulk is compact, region 
II disappears completely. In this case it is still important to consider particle production 
on the brane since otherwise the energy that was stored in the rolling tachyon might go 
into invisible closed string modes in the bulk, namely gravitons, which would not be an 
acceptable form of reheating. 



3. Gauge field solutions 



Our aim is to find out whether the energy stored in the homogeneous tachyon fluid can be 
efficiently converted into radiation, so that the universe at least has a long period of radi- 
ation domination before possibly giving way to the cold dark matter of the rolling tachyon 
condensate. We will do this by quantizing the gauge field in the tachyon background and 
computing the Bogoliubov coefficients that quantify the mismatch between the vacuum 
states of regions I and III (see fig. 1(b)); see for example That is, if we start in 

the vacuum state appropriate for region I, we find that it is no longer the vacuum in region 
III, and therefore radiation must be produced. 

The first step is to find the action for the gauge fields to quadratic order in the fields. 
Expanding ( ^.4| ) in the tachyon background described in the previous section, we obtain 



I dtdxd'^y < 



2 



f {dtAf - {VyAf - {d^Af, region I 

(dtA)^ - e-2*/<^ ({VyAy + {d.^A)A , region II 



yf^g-<7|x|/a /(5^^)2 _ (v^^)2 _ (5^,^)2 j ^ ^^^^^^ jjl 



(3.1) 

where yi are the coordinates of the large 3 dimensions, and we have absorbed the volume of 
any other compact dimensions which are merely spectators into the brane tension T. We 
employed radiation gauge {Aq = Vy ■ A = 0) and projected out the extra polarization by 
setting Ax = 0, which is consistent since this state turns out to have a mass gap in region 
III, unlike the massless components among the large three dimensions, A. (The factor 
g-2t/a -j-^ region II should really be (cosh(t/o))~^, but this not an important difference, in 
the spirit of the other approximations we have made.) In the following we will drop the 
polarization indices of the gauge field and write simply A instead of A. 

To make our analysis more tractable, we have ignored the time dependence of the 
metric due to the expansion of the universe in the above action. This neglect can be justified 
if the initial fast-roll regime of the tachyon, during which most of the particle production 
occurs, does not take more than approximately one Hubble time. We expect that this will 
be true if the string scale is somewhat below the Planck scale and strings are weakly coupled, 
since then H ~ y/2LT/M^ ~ 5,Mj(27rMp) [using the relation g^M^ = M^V/niT'f @ 
and eq. ( [4. 61) ], whereas the time scale for the tachyon roll is of order the string scale. 

In the previous section we gave detailed motivations for our choice of the ansatz for the 
classical tachyon background. It is worth emphasizing one further criterion: by choosing 
T{t, x) to depend only on t or a; in each region, we insure that the gauge field equations of 
motion can be solved using separation of variables, which greatly simplifies the task. 

3.1 Solutions in each region 

The solutions in region I are trivial since the background tachyon configuration is simply 
T = 0: 

A, = V \{a^e-'^-' + ale'^^') cos(A;^x) + (a^e"^"'"* + cile'^^') sin(fc^x) 

(3.2 



where km = rmr/L, uJ^ = p'^ + k^, and p is the momentum in the three large dimen- 
sions. We have spht the modes according to their parity in the extra dimension for later 
convenience. 

In region II, the equation of motion for A is 

A + u;^e-2*/M = 0, (3.3) 



which has the solutions 
1 



■t/a- 



COS 



1.x 



+ same with br, 



Cm, COs{kmX) Sm{kmX) 



(3.4) 



Near t = 0, these oscillate just like the region I solutions, but at large t the oscillations 
freeze and the solutions grow linearly with time. 

Region III is the important one at late times, since this is where the descendant brane 
and the standard model are supposed to reside. Here the equation of motion is 



(1 + g2) i^-A + V^Aj + A" - %gn{x)A' = 

using primes to denote dx- The solutions can be written as 
1 



(3.5) 



Ar 



E 



{dne-'^-' + 4e^^"*)/„(x) + (4e-'^"* + dle'^-') fn{x) (3.6) 



with 



UJr, 



fn{x) 



fn{x) 



No, 

N^f.q\x\/2a fcos{knX) 



1/2 



2kna 



sm 



^(l_e-.Wa)-l/2(l+g2)-l/4^ 



2a 
4L 



/4L 



1 + 



-1/2 



q\x\/2a sin(A;„x), 



(1+a-^^ 



n = 

n > 1 

n = 

n > 1 

n = 

n > 1 

n > 1 



(3.7) 



and kn = mxjL. 

These solutions have the desirable property, from the point of view of string theory, that 
there is a zero mode accompanied by a tower of heavy states [21, 24]. This is qualitatively 
similar to the spectrum of excited states of the open string, though we have sacrificed 
some of the similarity by taking the tachyon potential e~\'^\l°- instead of e~"^^/"^. With the 
latter choice one gets a more realistic spectrum of the form u)^ ~ + n/a^, which has 
the correct n-dependence to match string theory. The disadvantage is that the solutions 
in region II cannot be found analytically since A + a;^e~^*^/"^ A = is not a standard 
differential equation. We have therefore given up some of the quantitative similarities with 
the real theory for the sake of being able to go as far as possible analytically. 



3.2 Matching at interfaces 



To complete our task, we must relate the solutions in each neighboring region to each other 
at the interfaces t = and t = q\x\. This will impose relations between am, dm and bm, Cm 
and between bm, Cm and dn, dn. At i = the procedure is straightforward; A and OqA must 
be continuous, leading to the relations 



bm\ _ T^a^f^ 

Cm I ~ 



-Yi - iYo -Yi + iYo 
Ji + i Jo Ji — iJo 



(3.8) 



where the Bessel functions are all evaluated at aiVm- We used the Wronskian of Jq and Yq 
to obtain { \i.8\) . Notice that there is no mixing between different m values at this point. 

Matching at the t = q\x\ interface is more difficult. In this case we also demand 
continuity of A, but the fact that the prefactor y^l — (5^T)2 in the action is discontinuous 
means that derivatives of A are discontinuous as well. By integrating the partial differential 
equation for A along an infinitesimal path which crosses the interface perpendicularly, we 
can show that the following linear combination of derivatives must be continuous at the 
interface: 



A{FtA + qF^A') = from II to III 



(3.9) 



where Ft is the coefficient of A and F^ is that of A in the action. Namely, = 1 in 
region II and y^l + g2e-<?kl/a m, while F^ = e'^*/" in II and g-^l^l/^/y/l + g2 HI. 
This gives the rather cumbersome matching conditions 



= (6„ Jo(aw^e-'?l^l/") + CmYoiau;me-''\''\h) cos(fe„x) 

= (d„e-^^"^l"l + dte^^-"!^!) /„(x) (3.10) 



for A itself and 



4L 4::^ 



[bmJliaUJmC ^''''/") + CmYxiaUJmC '?l''l/")j COs(/CmX) 

- gfeme-"!^!/" f6„ Jo(aa;™e-^l^l/'^) + c^yo(au;™e-'?l^l/'^) ) sin(A:„x) 



+ 



^/TT^ f d^e-*^""!"! - die 



iu)„q\x\ 



fnix) 



(3.11) 



for the derivatives of A. Similar conditions hold among the odd parity modes [a, d ^ d,d, 
fn^fn, cos{kmx) sm{kmx)], which do not mix with those of even parity. 



3.3 Solution of matching conditions 

The technical challenge is now to solve for dn, dn in terms of bn,Cn- Normally this would be 



done by taking the inner product of each equation ( [3.10 , 3.11 ) with some function which is 



orthogonal all but one of the functions multiplying dn or di. But because of the diagonal 
nature of the II/III interface, the latter functions are products of orthogonal functions, 
which are no longer orthogonal. This makes it impossible to solve for dn , dn analytically. 

We should therefore choose some set of basis functions gi{x) and integrate them against 
the matching conditions ( 3.10| - 3Tl| ) to transform the latter into the discrete form 



EA.(y^i:«.™(t) (3.-) 

We were not able to identify any set of gi{x) that makes the matrix Lin even approximately 
diagonal (which would facilitate an analytic approximation), making a numerical solution 
necessary. After some experimentation, one realizes that the most efficient way to do this is 
to discretize the system on a spatial lattice at positions Xi, so that gi{x) = 6{x — Xi). This 
allows us to compute the matrices Lin and Rim without having to perform any integrals. 
An ultraviolet cutoff is thus introduced. Let i = 0, . . . ,iV so that Xi = iL/N. It makes 
sense to let the mode number of the spatial eigenfunctions fn{x) also range from to N; 
then ( 3.12| ) gives exactly as many equations as unknowns. The system can be solved by 



numerically inverting the matrix Lin [27|. The results are presented in the next section. 
4. Particle Production 

Our solutions for the gauge field in regions I and III are normalized so that am, aln, dm, dm 
are the correct creation and annihilation operators for particles in the distant past and 
future once the gauge field is quantized. We assume the universe starts in the vacuum 
state a„|0) = 0. But this state is not annihilated by d„, since the latter is a superposition 
of a, determined by the Bogoliubov coefficients a and /3, 



dn 
di 




(4.1) 



Therefore observers in the future will see a spectrum of particles in the final state given by 

Mn = {0\didn\0) = Y,\(^nm\^ (4-2) 



m 



In this section the numerical results for J\fn and the total energy density of produced 
radiation will be presented. 

4.1 UV sensitivity 

Ideally one should obtain convergent results in the limit that N ^ oo, where we recall 
that N is the number of sites of the spatial lattice in the extra dimension, introduced in 
the previous section. However we do not observe this from our numerical results; rather 



there is a steady growth with N. Our conjecture is that this is related to the discontinuity 
in the action which arises from the sign change in (d^T)'^ across the II/III interface. In a 
simpler situation where particles are produced due to a (spatially constant) time-varying 
background, the spectrum is proportional to the square of the Fourier transform of the 
background. If the latter has sharp features, the spectrum falls only as a power of energy, 
whereas a smooth background leads to exponential suppression of high energies. For exam- 
ple, the Fourier transform of e~*^ times a step function of time is (A -|- iu;)~^ whereas the 
Fourier transform of (l-|-t^A^) is proportional to e~^^^^^. In the former case, summing over 
high frequency modes could lead to nonconvergent results. We expect the high frequency 
contributions to be suppressed by a factor like e"''^'/^ if we had a more realistic tachyon 
profile whose derivative changed smoothly over a time 1/A. We have not yet found a way 
of altering T{t, x) to incorporate this behavior while still allowing us to solve for A{t, x) 
analytically. 

To remove sensitivity to the lattice spacing, we therefore try to model the expected 
effects of smoothing out the background tachyon solution by inserting a convergence factor 
by hand, so that ( |3.12| ) becomes 

Once this is done, the limit ^ cxd is well-behaved. We have thus essentially traded 
the original ultraviolet cutoff N/L for a new one. A, whose physical meaning is more 
transparent. Figure 2(a) shows the dependence on of the low-momentum zero- mode 
production, J\fo{po), at po = 0.05/a, for several values of A and q (recall that the latter 
parameter determines the slope of the interface between spacetime regions II and III). The 
convergence with A^ is faster for smaller values of q; numerical limitations therefore prevent 
us from accurately studying large values of q. Fig. 2(b) shows the dependence of zero- 
mode production on moderate values of q; large values of q at fixed A^ give exponentially 
increasing results due to the term Yi{auJrnS'~'^^^^) in ( |3.11 ), but these nevertheless become 



well-behaved again as A^ is increased to sufficiently large values. On physical grounds one 
might anyway expect q < 1 due to the finite speed of propagation of the kink. 

4.2 Spectra 

We note that each state in region III is distinguished not only by its mode number n, but 
also its momentum p in the large dimensions, which we have until now suppressed except 
for its appearance in the energies tD^, eq. ( |3.7D . In figure 3 we show the dependence of 
■!^n{p) on these two variables. The spectrum of zero modes falls monotonically with p as 
expected. Although the nonzero mode spectra temporarily rise with p, their contributions 
to the total energy density are smaller than that of the zero mode. We have taken ^ = 1 
and A = 1/a for the parameters controlling the steepness the kink in spacetime and the 
suddenness with which it forms. The dependence on the size of the bulk L can be seen in 
fig. 3(b). 

Here and in the remainder we have made the simplifying approximation of ignoring 
the odd-parity modes and keeping only the even ones, which include the dominant zero 



mode. We expect that this gives a shght underestimate of the actual efficiency of particle 
production, by no more than a factor of order unity. 




0.4 



0.2 



N= 100 
N = 200 



Figure 2. (a) Dependence of A/odJo) on forpo = 0.05/a, L = 2a, A = 0.5/a, 1/a and (7 = 0.5, 1, 2. 
b) Dependence of Afo{po) on q for N = 100, 200, at po = 0.01/a, L = 2a and A = 1/a. Physically 
meaningful results correspond to the N 00 limit. 



To find the total energy density of produced radiation, we should sum over both n and 



p: 



dpr 

dp- 
dp 



d^p 
(2^ 



(4.4) 



The heavier modes are counted because they will presumably decay very quickly into 
massless standard model particles. In figure 4 we show the differential energy density, 
as a function of momentum, for several values of the other parameters. 
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Figure 3. (a) Momentum dependence of the spectrum of produced particles, Mn{p) for the string 
excitation quantum numbers n = 0, . . . , 5. (b) Afn{p) as a function of n for L — 8a, 4a, 2a, a and 
p = Pq = Q.Ol/a. Both graphs are for q = 1, A = 1/a . 
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4.3 Energy density produced 

We turn now to the main results, the total energy density pr of produced radiation, ob- 
tained from integrating eq. ( |4.4| ). In figures 5(a-c) we graph the dimensionless combination 
ln(a'^/3r-) as a function of the main unknown parameters, g, A and L. It is clearly an in- 
creasing function of all three parameters. The "critical value" shown in these figures, which 
we would like pr to exceed, is derived in the next subsection. We will argue that this is 
the value at which reheating starts to significantly deplete the energy stored in the rolling 
tachyon fluid. 




Figure 5. (a) ln(a'*pr) versus aA for L = a, 4a, 8a. (b) \n.{a'^pr) versus q for several values of L and A. 
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4.4 Efficiency of reheating 

Our goal is to see if the produced pr is a large enough fraction of the energy density which 
is available from the homogeneous rolling tachyon fluid, Pt- We thus need to specify the 
value of Pt which is predicted by string theory. In terms of the effective tension of the 
original brane and the size of the extra dimension, 

Pt = 2LT (4.5) 

By the effective tension, T, we mean the energy density in the (4+l)-dimensional spacetime 
which includes x. If we started from a nonBPS Dp-brane, then T = ^plT^y where V is 
the volume of compact dimensions within the brane excluding x, and 7^ is the tension of 
a BPS Dp-brane, 

r = l(2^)-P(V^)-(P+i) (4.6) 

Therefore pr depends on gg, V, L and a' . However not all of these are independent; they 
are related to the fine structure constant of the gauge coupling evaluated at the string scale 



= 2i2LV) ^'-'^ 



Interestingly, if we use this to eliminate the V dependence from ( [4.5D we find that the 
dependence on the string coupling and on L is also gone: 

Pt = r- (4.8) 

^ v^(2^)3a(M,) ^ ' 

To find out how much energy is available for reheating, we must also consider the 
descendant brane's 4-D energy density, p^. This is obtained by integrating over the extra 
dimension in ( |2.4| ) for the kink profile T = qx and taking the limit q ^ oo. The result is 
simply Pf = 2aT . Therefore the excess energy density which can be used for reheating is 

/\p = p^-pf = p^ [l - I) (4.9) 

Hence we should demand that L > a to get any reheating at all. 



In order to compare to the energy density of produced radiation, we need to 
know the parameter a in terms of the string length. In our simphfied model of tachyon 
condensation, this can be determined by demanding that the tension of the descendant 
{p— l)-brane match the string theoretic value for a BPS state, given by (4.6) with p — > p—1. 
In the limit (7 — > oo, the ratio of the initial and final brane tensions in our field theory 
model is (limg^oo \/l + (7^ e~'l^'l/"dx)~^ = l/2a, whereas the string theoretic value is 
\/2Ms/(27r). We therefore obtain 

Now we can quantify the efficiency of reheating. Our results for are expressed in 
units of a~^. Using ( [4.10| ) we can convert the available energy density Ap into the same 
units to find the critical value of pr, call it pc, for which the conversion into radiation would 
be 100% efficient: 

Pc = ^ a-^ ^ l.la-^ (4.11) 

32V2 a{Ms) 

We take the fine structure constant to be 1/25, and we omit the factor (1 — a/L) since it 
would require fine tuning of L to make it very relevant. The figure of merit for reheating 
is therefore Pr/pc- The critical value pc appears as a dotted line in figures 5(a-c). We see 
that for large enough values of A, q, or L, reheating can be efficient. 



5. Conclusion 

We have made a case for reheating the universe after brane-antibrane inflation by pro- 
duction of massless gauge particles, due to their coupling to the fast-rolling tachyon field 
which describes the instability of the initial state. Our results are encouraging, indicating 
that if the extra dimensions within the original branes but transverse to the final one are 
large enough compared to the string length scale V^, reheating can be efficient enough 
so that radiation dominates over cold dark tachyon matter. Depending on details of the 
tachyon background and the compactification, an extra dimension of size ^ 10\/^ could 
be sufficient. (For example, with L = 4a in an orbifolded compactification, the size of the 
extra dimension would be ATT/y/2M~^ = 9Vo/. In the present analysis we counted only a 
single polarization of one U{1) gauge boson and ignored the odd parity modes in the extra 
dimension, so the real situation could be less constrained. 

The calculation is complicated, and we have made it tractable by invoking a number 
of simplifying assumptions. We considered formation of a kink in the tachyon field rather 
than a vortex. We used a simplified version of the tachyon action rather than the one 
which has been derived in BSFT. We used an ansatz for the background tachyon field 
which is a good approximation to the actual solutions in the regions where it depends only 
on space (the kink) or on time (the homogeneous roll), but we do not know how suddenly 
it makes the transition between these two regimes, hence we parametrized this uncertainty 
by introducing a cutoff A on the tachyon field's derivative. We have also left the slope q 
of the tachyon kink profile T = qx, during kink formation, as a free parameter. Moreover 
we have ignored the expansion of the universe, which is only correct if the brunt of the 



tachyon roll completes in a time not much exceeding the Hubble time. We have also ignored 
the back-reaction of the produced particles on the tachyon background, so our criterion of 
requiring the produced energy density to meet or exceed that which is initially available is 
a crude one. 

To do a better job, more numerical computations will probably be necessary since it is 
hard to solve the gauge field equations of motion in a tachyon background that depends on 
both space and time. One check that might be carried out relatively easily is to numerically 
obtain the time- and space-dependent solution T{x,t) for the more realistic action and 
compare its features to those we have assumed. This is in progress. Very recent work on 
how to generate time-dependent tachyon solutions in ref. [^] might prove to be helpful 
here. 



As we were finishing this paper, a related work appeared |29]. The latter assumes that 
there is a constant rate of decay of the tachyon fluid into radiation. We have presented 
arguments to the contrary. If such an assumption could be justified, it would greatly 
enhance the viability of inflation via brane-antibrane annihilation. 

We thank C. Burgess for useful suggestions during the course of this work. We also 
thank D. Easson, K. Bardakci, A. Sen, P. Kraus, P. Steinhardt, and H. Tye for helpful 
comments and correspondence. 
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